Convergence of Perturbations for a Big Bounce in Loop Quantum Cosmology 
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We investigate the convergence behaviors of the scalar and the vector perturbations for a big 
bounce phase in loop quantum cosmology. Two models are discussed: one is the universe filled by a 
massless scalar field; the other is a toy model which is radiation-dominated in the asymptotic past 
and future. We find that the behaviors of the Bardeen potential of the scalar mode near both the 
bounce point and the transition point of the null energy condition are good, moreover, the unlimited 
growth of the vector perturbation can be avoided in our bounce model. This is different from the 
bounce models in pure general relativity. And we also find that the maximum of an observable 
vector mode is inversely proportional to the square of the minimum scalar factor Obounce- This 
conclusion is independent with the bounce model, and we may conclude that the bounce in loop 
quantum cosmology is reasonable. 

PACS numbers: 98.80.-k,98.80.Cq,98.80.Qc 



I. INTRODUCTION 

There are several ways to solve the singularity problem 
in cosmology one of them is the bounce model 
In the studying of the bounce model in pure general rel- 
ativity (GR) 5], the scalar hydrodynamic perturbation 
can lead to a singular behavior of the Bardeen potential. 
The further study show that the divergence of the 
Bardeen potential can be avoided only in some models 
with special matter. 

Another mode of perturbations must be considered is 
the vector perturbation. It is known that the vector mode 
will decay quickly in expanding phase of the universe. So 
it exhibits a growing mode solution in a contracting uni- 
verse foj . This growing will in general lead to the break- 
down of perturbation theory near the bounce. So, it is 
necessary to check the behaviors of the vector perturba- 
tions near the bounce. 

Generally speaking, the bouncing phase originates 
from a quantum effect of gravity. So it is interesting 
to study the bounce model in quantum gravity theory. 

At present, the issue of finding a complete theory of 
quantum gravity is still open. In current approaches, 
one of the most active is l oop quantum gravity. Loop 
quantum gravity (LQG) |lC)| - [l2| is a mathematically well- 
defined, non-perturbative and background independent 
quantization of general relativity. And, Loop Quantum 
Cosmology (LQG) [I^, a symmetry reduction of LQG to 
the homogeneous and isotropic spacetime, has achieved 
many successes. A major success of LQG is the resolution 
of the Big Bang singularity [l^ - [T6j . this result depends 
crucially on the discreteness of the spacetime. Instead of 
Big Bang, there will be a Big Bounce. 

In GR bounce, the divergence of the Bardeen potential 



occurs at two point i.e., near the bounce point and 
near the transition point of null energy condition (NEG). 
The difference is that the LQG bounce is governed by a 
discrete quantum geometry [l7|. This will lead to some 
different behaviors of the Bardeen potential. In this pa- 
per, we consider the behaviors of the perturbations near 
both the bounce point and the transition point of NEG 
under the framework of the effective theory of LQG. 

The paper is organized as follows. In Sec. |lll we give 
the framework of the effective theory of LQG with holon- 
omy corrections, it can yield the bounce background. In 
Sec lIIIl we introduce the scalar perturbation based on 
the Sec. |lll Two models are analyzed in this section: 
one is the universe filled by a massless scalar field; the 
other is a toy model which is radiation-dominated in the 
asymptotic past and future. In Sec. |lV]we discuss the 
vector perturbation near the bounce. The discussion and 
conclusions are presented in Sec. |Vl 



II. BACKGROUND OF LQC BOUNCE 

The canonical variables used in LQG are the Ashtekar 
connection A^^ and the densitized triad Ef p^| - [T^ , where 
= -I- Kl with r the spin connection and K the ex- 
trinsic curvature, and Ef = e'^/\dete\\ with e"e^ = q°-^ 
and qab the spatial metric. For a spatially homoge- 
neous and isotropic universe model (FRW metric), the 
Ashtekar variables can be reduced to the diagonal form, 
i.e., Aj, = c5l and Ef = p5°: [H. Therefore, the basic 
canonical variables for the gravitational field are (c,p) 
and for the scalar field {(p^p^). Here we denote the back- 
ground variables with a bar. In this paper, we consider 
only the flat space universe. Thus the canonical variables 
can be expressed in terms of the standard FRW variables 
as: (c, \p\) = {-fa, a?), where a is the scale factor; and the 
effective Hamlitonian of the considered model is given by 
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where the factor 7 is called the Barbero-Immirzi param- 
eter which is a constant of the theory. 

In the process of quantization, we can find that there 
is no operator corresponding to the canonical variable 
c itself but we can return to the holonomy. This fact 
can lead to the so-called holonomy correction in an effec- 
tive theory of LQC. The effects of this correction can be 
obtain b y si mply replacing the c to sin(/2c)//2 with the 
choice of[2a| 



where the prime denotes the derivative with respect to 
the conformal time 77 and 7J? = c. In addition, the Klein- 
Gordon equation can also be derived as follows: 



sin(2/i7.ft) 



27/i 



(9) 



From Eq. (|4]) and the relation between cosmic time and 
conformal time, one can get the motion equation of p 
with conformal time 



(2) 



p' = 2p 



sin(2/i7.S) 



27/i 



(10) 



where A = 2\/3njlp is a area gap, and Ip denotes the 
Planck length. So, the effective Hamiltonian with holon- 
omy correction is 



H, 



eff 



sin(/ic) 



+ 



1 



Pi 



2 \p\V- 



.+\p\y^vm. 

(3) 



From now on, we focus on a positive p. 

The equations of motion can now be derived by the 
using of the Hamilton equation 



(4) 



So, one can further define a conformal Hubble parameter 
H by 



2p^ ' 



27A 



Moreover, we can also define 



©1 



n 2 



62 := H'-I'-p— 
op 



(11) 

(12) 
(13) 



where the dot denotes the derivative with respect to the 
cosmic time t and the Poisson bracket is defined as 



{f,9} = 



inGj 



3 

df dg 



dc dp dp dc 
df dg 



dip dp,f, dp,f, dip 



(5) 



From this definition we can obtain two elementary 
brackets 



871^7 



(6) 



By using these brackets, one can derive two equations, 
i.e., an effective Friedmann equation and an effective 
Raychaudhuri equation. However, in theory of pertur- 
bation, the use of the conformal time 77 may be conve- 
nient than the cosmic time. The conformal time can be 
related to the cosmic time t through the scale factor a, 
adr] = dt. Thence the effective Friedmann equation and 
the effective Raychaudhuri equation with conformal time 
are respectively [23^ . 
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In fact, 61 and 62 can be taken as two effective correc- 
tions to the evolution equations of the Bardeen potential. 
This is because that one can obtain the classical evolu- 
tion equations by taking Si — ^ and ©2 — >■ at the 
same time. 

We rewrite the Eqs.© and ([S]) in terms of 61, ©2 and 



M^-e, = ^pxps, (14) 
-i(H2-6i)-|(H'-62) = ^pxP„ (15) 



where 



and 



2p 



2p 



(16) 



(17) 



are the energy density and pressure of scalar field respec 
lively. 

By using Eqs.lfTTj) and (fT6|). one can rewrite the Eq.(fT4] 

to 



ir2 _ j2 



p p,[l 

Pc 



Ps 



where 



^272 A ■ 
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Obviously, one can easily check that when ps = pc, 
EI = 0, and this means that a bounce occurs. The bounce 
density pc is related to A which is the smallest eigenvalue 
of area operator, so this LQC bounce is originated from 
the quantum effects of spacetime. 

From Egs- lfM)) and (|15|) we can get a useful relation 
equation 



P. = 



1 



where 



61 + 62. 



(20) 



(21) 



Thus, the relationship [2l| between the null energy con- 
dition and the stress-energy tensor can be written as: 



NEC ^ 



(22) 



III. SCALAR PERTURBATION WITH 
HOLONOMY CORRECTIONS 



In this section, we introduce the scalar perturbation 
based on the Sec. [IT] and analyze in detail the follow- 
ing two models: one is the universe filled by a massless 
scalar field; the other is a toy model which is radiation- 
dominated in the asymptotic past and future. In our 
models, the spacetime is described by the metric 



= a'^irj) [-(1 + 2$)d772 + (1 - 2'^)Sabdx''dx''] , 

(23) 

where $ and ^ are the Bardeen potential. In the case 
of vanishing anisotropic stresses, $ and are equal (2^ . 
Therefore, we set $ = from now on. The evolution 
equations of the Bardeen potential with holonomy cor- 
rections have been given in [2^: 
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$" + i 3H + 4 
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2H2 + 4H' - m' - 3p^ 
op 



7M 



(f> = AnGp Sps 



(24) 



$ = iirGpSPs 



(25) 



By Using the definitions of Egs. p^ and p^ . we can 
rewrite the Eqs.(l24l) and ([25]) as 



dps 



dS 



Ps 



For the hydrodynamic matter, T can be interpreted as 
the sound velocity. In this paper, we focus on a adiabatic 
perturbation only, so we have 



T 



SJ^^P^ 

SPs P's 



(29) 



From Egs. p^ . (ITS|) and (PT|) we can obtain 



- 3H$' - [2tf + H' - 46i - 62] $ = 47rGp Sp^, 

(26) 

$" + ^3H + $' + [2tf + H' + 362] $ = 4:TrGp SP^. 

(27) 

We want to get the evolution equation of the Bardeen 
potential but the matter is not contained in it. So we 
should obtain the relationship between dps and SPs- In 
general, the pressure perturbation can be separated into 
two parts of the adiabatic and the entropic perturbation 
as follows 

= (^] Sps + (^] SS = TSps + tSS. (28) 



T 



1 / 2(3' + &[-2m&2 
3 V 2H/3 + G[ - 2H62 



(30) 



Using Eq.dini) and inserting Eq.(l27l) into Eq.^^, one 
gets 



= 



+ [(1 + T) (2H2 + H') - 4T6i + (3 - T) 62] 

(31) 

The equation for the mode of wave number k is 



= 



3H(1 + T) + |^ 



-f (1 + T) (2H2 + H') - 4T6i + (3 - T) 62] 

(32) 

A. Massless scalar field 

In this subsection, we discuss a simple model which the 
universe is filled with a massless scalar field, i.e., V{ip) = 



4 



0. The exact solution oi p(t) is [2J| 
where 



ill 



8^3 



-,1/3 



is the minimal value of p(t), and 



(33) 



(34) 



(35) 



A = 127rGp;. 

Because V{ip) — 0, one can get p^p — 0, p^p — const., then 
Pmin and ^ are constant too. 

In our discussion, we need to know the evolution func- 
tion of p with rj. It can be obtained by using the relation 
of adrj = dt. So, we have 



■n 



1 



p(t)l/2 



dt. 



(36) 



However, the result of the integration is complex. To 
keep our discussion simple, we consider a asymptotic be- 
havior of p{t): 



A'^h^^\ \t\ 



oo, 



Itl^O. 



{far from the bounce) 
(near the bounce) 
(37) 

From the integration of Eq. p6|) . we can get the asymp- 
totic behavior of 77: 



= <^ 2 



A 1/6^2/3^ _^ (ij^^ from the bounce) 



0. 



{near the bounce) 
(38) 

Then, we have an approximate relation between 77 and t 



(39) 



Figdlshow that \t]\ = {'S/2)A^'^/^t'^^^ is a good approx- 
imation for 77 = J{l/a)dt. The most straightforward way 
is to insert this approximate relation to Eg. ([55]) . How- 
ever, this relation is only a asymptotic behavior corre- 
sponding to the case of |t| — > 00, not the behavior on all 
time. So we should consider the asymptotic behavior of 

p{v)- 

Under the approximation of Eq. p9p . the asymptotic 
behavior of 77(77) is 



t\ — >■ 00, {far from the bounce) 
Itl^O. 



{near the bounce) 
(40) 

Now, we construct a function of 73(77) approximately. 
The approximate function should satisfy the asymptotic 
behavior of 77(77) . One class of such functions are 



p{v) 



1/n 



(41) 



a 

O 

c 
o 




Cosmic time 



FIG. 1: The relation between the conformal time ri and the 
cosmic time t. The solid line is 77 = j{l/a)dt, and the dash 
line is 7? = (3/2)^1-^/^*2/3^ 



where 77 is a natural number. 

One can find that if n is even, there will not be an 
absolute value like I77I appeared in the equation. So, from 
now on we set 77, = 2. There is another reason to choose 
77 = 2 that only 71 = 2 can lead to a evolution equation 
of the Bardeen potential which have analytical solution 
(see Eq.dUl)). 

So we set the approximate function of 77(77) is 



p{v) 



[Pn 



Arj 



2\l/2 



(42) 



where A = ^A. Eq.(|31]) is different with Eq.(|331)- Equa- 
tion (|33p is exact evolution of 77 but Eq. p2l) is a approx- 
imation of Eq. ([551) . 

Under the condition of V{ifi) ~ 0, we have Pg — Ps, so 
T = 1. From Eqs.dH]), ^ and ^ we can obtain 

/? = 3(h2-6i), (43) 
62 = 2H2 - 26i +H'. (44) 

Using these equations and Eq. ([5^ one can get 



(P + 46i)$fc =0. (45) 



Using Eq. (|42|) and Eq. (|45l) . we can discuss the behaviors 
of the evolution equation of the Bardeen potential in the 
following cases. 



1. Near the bounce 

Near the bounce point, |7;| — ^ 0, and the leading order 
of the coefficients of evolution equation are 



8H- 



.61 



1 



46 1 



(46) 
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And the evolution equation changes to 

and the solutions of Eq. (|47l) is 

<fk{v) = V^iCMICv) + C2Z_,(/Cry)] (48) 

where 

and Zi, is Bessel function, Ci, C2 are arbitrary constants. 
The leading order of this solution is 

*fc(r/)«$(l)+$(2)?7'^ (50) 

where <i>(i) and $(2) are constants which related to k. 

One can found that, v < 0, and the second term in 
Eq.([5(Il) is divergence. But we can choose the arbitrary 
constant C2 — 0, which means $(2) = 0. Thus we can ob- 
tain a convergence solution of scalar perturbations near 
the bounce. 



In GR, if the matter always satisfies NEC, there will be 
no bounce. The reason is that the GR bounce is led by 
some exotic matter which violates the NEC. However, in 
LQC, the bounce is originated in the discrete spacetime 
geometry. It is a quantum effect of spacetime. So, even if 
the matter never violates the NEC, there is also a bounce 
phase. 



B. Toy model 

In this subsection, we extend our discussion slightly, 
and consider a toy model which was introduced in |18l . 
[Tgj . In this model, the universe is radiation-dominated 
in the asymptotic past and future. In other words, the 
asymptotic behavior of 0(77) should be a oc 77 or p cx 77^. 
So we assume that there are some V^ip) that can make 
the form of ^(77) as[29| 

p(^) = (6'"-^a7;2™)i/™, (51) 

where b, a are constants, and & > 0, 771 > 1. 

1. Near the bounce 



2. NEC transition problem 

In our model, the point of NEC transition is obtained 
form j3 = 0. When V{^) = 0, the /3 is Eq.gSl). From 




-3 -2 



-1 1 
Conformal time 



FIG. 2: The evolutions of 
point. 



-61 and l3 near the bounce 



the Figl2l we can see that (3 is always positive in our 
model. So there is no NEC transition. And then, there 
is not the problem of the divergence near the point of 
NEC transition. 



Under the assumption of Ea. ((5T|) . we can obtain the 
evolution equation of the perturbation near the bounce 
point from Eq. ((5^ 



2m -2 , ^, r- 

1 (2777- l)(V6 + 5) 



3a 



In fact. Eg. ([5^ is only accurate to its leading order. We 
can obtain the following analytical solution 

M^l) = 77(l-™)''[$(i)Z.(^7/-") + $(2)Z_,(^7/-™)], 

(53) 

where is the Bessel function, ^j-i) and $(2) are con- 
stants, and 



1 



1 — 77^ 
ik 



777,-1 



'6(2™"l)/2(2771-2)(l + \/6) 



777 — 1 y 3a 
The leading order of Eq. (l53l) is 

$fe(77)«^"^[$(i)+ $(2)772(1-")-] 



, (54) 



(55) 



(56) 



Because 2{l-m)v = {2m-2)/^/b+{m-l){^/b + h) > 0, 
so the behavior of the perturbation near the bounce point 
is good. 
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2. NEC transition problem 

Near the bounce point, the leading order of /3 is 
^ a(2m-l)(Vb + b) ^ 



(57) 



So, there must be a time point corresponding to /? = 0. 
To get this point, we consider the next order of /3 

_ a(2m-l)(Vb + &) 
, a%2m - l)(4m - + b) + b] 



(58) 



Near the bounce point, the leading order term is rj^"^~'^ 
term, so /3 < 0. When jryj get larger, the leading order 
term will be 7y4m-2^ so /? > 0. The transition point is 
/3 = 0, it leads to a transition point 



± 



a(4m - 1) 



(59) 



770 denotes the transition time. 

When we discuss the perturbation near the t^o, we can 
shift the origin of the time such that 770 = 0. If we do 
that, the function of p{vi) will be 



p(r;) = [6" + «(r; + r;o)^"]^/". 



(60) 



Near the point 770, ~ and T « —2/3. Moreover, from 
the definition of /3 and /3 = 0, we have ©2 = ei-H^+H'. 
Thus, Eq.(1211) changes to 



61 



■r 



56i - 



= 0. 
(61) 



Near the ryo, the leading order of this equation is 

+ (^-^fc2 + L2^ $fe = 0, (62) 

where Li and L2 are constants which can be related to 
TO, A etc. The solution of this equation is 



*/c = ^'(1) exp 



-(Li + k)r] 



+$(2) exp 



1 



(Li - k)r] 



where ^{i) and $(2) are constants again and 




L, 



Ll 



(63) 



(64) 



Near the NEC transition point, we can see a convergence 
behavior of the perturbation. 

IV. VECTOR PERTURBATION WITH 
HOLONOMY CORRECTIONS 

In the most of research on perturbations of cosmology, 
a lot of attention have been focused on the scalar and 
the tensor mode. The reason is that the vector mode 
will decay quickly in expanding phase of universe. 

However, in the pre-bounce phase of the bounce mod- 
els, the universe undergoes a contracting. It is shown 
that, in contrast with the expanding phase, the vector 
mode will exhibit a growing behavior [9| . The unlimited 
growth of the perturbation will breakdown the perturba- 
tion theory. So, it is necessary to check the behaviors of 
the vector mode near the bounce. 

The effective linearized equations of vector mode with 
holonomy corrections have been given in [25| (in Newton 
gauge): 



2a 



^V25^ = (87rG)2(p + P)V\ 



27/iJ? 



(65) 
(66) 



where 5* and V* are the metric perturbation and the 4- 
velocity perturbation of the vector mode respectively, 11* 
is the anisotropic stress, and C/j"' is the anomaly term 
[2^. To have a consistent set of the evolution equations, 
we require the anomaly term to vanish i.e. G^^' — 0. 



S in Fourier mode of k: 



VI 



2(8^G)2p(p + P) 



(67) 



From Eq. (j65p . one can obtain a relation between V and 



If we do not take into account the anisotropic of per- 
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turbations, the Eg. (1661) will be: 



1 



Si 

2a 



1 



sin(2/i7^) 



and the equation of the Fourier mode k is 

sin(2/i7j|) 



dr] a 



SI = 0. 



= 0, 
(68) 

(69) 



We will solve the Eq. (pU|) in two models which is the same 
as the Sec. |llTl 



Massless scalar field 



This means that the limited growing is the same as 
Eq. ([7T|) . and we also have 



(6™ + Q;r;2™)i/™ ■ 



(75) 



We find that, the maximum of (p + P)Vl is also in- 
versely proportional with b ~ O'lounce- 1^ means that the 
maximum of the observable quantity {p + P)Vl. near the 
bounce is inversely proportional to the square of scale 
factor at the bounce point, and this conclusion is inde- 
pendent with the model. 



For the same massless scalar field model discussed in 
Sec. IIII A[ the evolution of the scale factor is Eg.ip^. 
and the Ea.(|Ml) changes into 



r) 

^Si + ^r^Sl = 0. 

OTJ Pmin 

We only consider the leading term and thus 

A2 



S]. (X exp 



-■n 



(70) 



(71) 



We can see that, even though the vector mode is grow- 
ing when the universe contracting to bounce, there is a 
maximum at the point of the bounce. It means that the 
vector mode should not growing unlimited. 

As pointed out in |9|] that, only the combination 
(/9+P)V* appears in the energy momentum tensor; there- 
fore it is this combination that could in principle be ob- 
servable and may thus be called physically relevant. 

From Eq. (|67|) . we can obtain: 



exp 



2p„ 



(72) 



This also have a maximum at the point of bounce and it 
inversely proportional with pmin = aL«nce- 



B. Toy model 



For the same toy model discussed in Sec. IIII Bl the 
form of p{r]) is taken as Eq. ((ST|) and the Eq. (IMl) approx- 
imating into the leading order becomes 



2a 
6™ 



We can obtain 



SI oc exp 



2m 



(73) 



(74) 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we examined the behaviors of the scalar 
and the vector perturbations in the bounce phase of the 
effective theory of LQC. Differing from the bounce model 
in [21| , the scalar perturbations in our model is not diver- 



gence near both the bounce point and the NEC transition 
point. Another conclusion is that the vector mode of per- 
turbations have maximum at the bounce point, and this 
maximum is inversely proportional to the square scale 
factor at the bounce point. 

In the model of GR bounce, the emergence of bounce 
phase is rooted in the matter in the universe. According 
to the singularity theorems [1^ [l^] , if one requires the 
matter satisfies the energy conditions, the universe can 
be emerged from an initial singularity. However, there 
is no evidence that the exotic matter which violates the 
NEC does not exist. So one can choose some exotic mat- 
ter to make the universe to experience a bounce. It is 
also because of this, the behavior of the bounce and the 
perturbation near the bounce point is decided by some 
exotic matter which have been chosen. Therefore, we can 
select the matter carefully to make the behavior of per- 
turbation near the bounce have good performance like 
in 0-11] • But too much artificial factors will make the 
physics of the model unnatural. 

On the other hand, the LQC bounce is originated in 
the discrete spacetime geometry. Just like the model in 
Sec. IIII Al even if the matter satisfies the NEC, there 
is also a bounce phase. So the behavior of the bounce 
is decided by the effects of discrete spacetime geometry. 
From the analysis in Sec. IIII Al and IIII Bl one can find 
that, the effects of discrete spacetime geometry lead to 
the convergence of the Bardeen potential. 

One should note that, our discussion is in the frame- 
work of the effective theory of LQC, so we find that this 
effective theory reflects the nature of quantum spacetime 
geometry effectively. Moreover, from the discussion of 
this paper, we also can obtain the conclusion that the 
existence of LQC bounce is reasonable, it do not lead to 
unbounded growth of the perturbation. 
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